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Abstract. An addition and product formula for the Hahn-Exton g-Bessel function, previ- 
ously obtained by use of a quantum group theoretic interpretation, are proved analytically. 
A (formal) limit transition to the Graf addition formula and corresponding product formula 
for the Bessel function is given. 



1. Introduction 
A classical result for the Bessel function 
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^ ■ \ye % y\ < \x\, due to Graf (1893), cf. [19, §11.3(1)], for general u, and due to Neumann 

(1867) for v = 0, cf. [19, §11.2(1)]. In case v G Z the conditions on x, y in (1.2) can be 
removed. 

Several g-analogues of the Bessel function (1.1) have been studied. The oldest q- 
analogues have been introduced by Jackson in 1903-1905, cf. the references in Ismail 
[4]. The g-Bessel function studied in this note is the so-called Hahn-Exton g-Bessel func- 
tion which has been introduced by Hahn (1949) for a special case and by Exton (1978) in 
full generality, cf. references in Koornwinder and Swarttouw [13]. 
The Hahn-Exton g-Bessel function is defined by 

(1-3) Uz- q) = z (g;g)oo KPi { qa+1 ; q, QZ 
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Here q G (0,1), (a; q) = 1, (a;g) fe = rT*=o C 1 ~ a Q % ) for k E~N, (a; g)^ = lim^^a; q) k 
and the g-hypergeometric function is denned by 

/ai, a 2 , . . . , a r \ = y, (ai; g) fc (a 2 ; g) fc . . . (a r ; g) fc //^ifc^-iA ( s " r+1 ) 
V &i, / ^ {q;q)k{bi;q)k---{b s ;q)k V / 

The notation for g-shifted factorials and g-hypergeometric series is taken from the book 
by Gasper and Rahman [3] to which the reader is referred for more information on this 
subject. 

The goal of this note is to prove the formula 
(1 4) 00 

1 ' } = ^(^^ (x+y+fc) ;?)J,+ fc ( J R^ (y+fc+iy) ;?)J,(^ (2 - fc) ;?). 

k= — oo 

This formula is valid for z e Z, R, x, y, v G C satisfying \r\* < l, ft( x ) > _1, 

and R ^ 0. 

The formula (1.4) has originally been derived for u, x, y G Z, R > by Koelink using 
the interpretation of the Hahn-Exton g-Bessel function as matrix elements of irreducible 
unitary representations of the quantum group of plane motions, cf. [7, §3.5]. This quantum 
group theoretic interpretation of the Hahn-Exton g-Bessel function is due to Vaksman and 
Korogodskh [16]. The quantum group theoretic derivation of (1.4) in [7] is modelled on 
the group theoretic derivation of Graf's addition formula (1.2) as presented by Vilenkin 
and Klimyk [18, §4.1.4(2)], so we call (1.4) a g-analogue of Graf's addition formula for the 
Hahn-Exton g-Bessel function. 

In section 4 we present a (formal) limit transition of (1.4) to (1.2) as g tends to 1. The 
approach employed is based on a theorem presented by Van Assche and Koornwinder [17] 
which has been used to show that the addition formula for the little g-Legendre polynomial 
[12] tends to the addition formula for the Legendre polynomial. The case v = 0, R = 
of the addition formula (1.4) can be obtained by taking a (formal) limit in Koornwinder 's 
[12] addition formula for the little g-Legendre polynomials using the limit transition of the 
little g-Jacobi polynomials to the Hahn-Exton g-Bessel function, cf. [13, prop. A.l]. 

Kalnins, Miller and Mukherjee [6] have given another derivation of the addition formula 
(1.4). They consider representations of the Lie algebra of the group of orientation and 
distance preserving motions of the plane. Instead of exponentiating the representations 
using the exponential function, they use a g-analogue of the exponential function. In 
a particular case the matrix elements can be expressed in terms of the Hahn-Exton q- 
Bessel functions. They give a decomposition of the tensor product (not the standard 
tensor product, but one closely related to quantum groups) of two representations and 
the corresponding Clebsch-Gordan coefficients are in terms of the Hahn-Exton g-Bessel 
functions. Comparison of matrix coefficients yields the addition formula (1.4). 

More results on Graf type addition formulas for the Hahn-Exton g-Bessel function, the 
Jackson g-Bessel function and other g-analogues of the Bessel function can be found in [1, 
2, 5, 6, 7, 8, 9, 13]. 

We start the proof of the g-analogue of Graf's addition formula (1.4) for the Hahn-Exton 
g-Bessel function by proving the product formula 

(-l)-g-l- J m (Rq^+y); q)J„_ rn (Rq?(y+»- m >; g) 
(I 5) 00 

= J2 q z Uq h{m+z) ,q)J.-M h ,q)URq Uv+v+z) ,q) 
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valid for m G Z, R, x, y, v G C satisfying > -1, g i+»0«0+»(v)|.R|2 < 1 and # ^ 0. 

This is done in §2. The proof of the addition formula (1.4) is a direct consequence of (1.5) 
as is shown in §3. The product formula in the case v = m = 0, R = q~? is mentioned 
(without proof) by Vaksman and Korogodskh [16, p. 177]. 

In §5 we will show that (1.5) is a (/-analogue of the product formula for Bessel functions, 

(1.6) J v+m (x) J m (y) = ± J v (aA 2 + V 2 ~ Ixy cos ^) ( *~ ^ ) ' e"^ #. 

The product formula is a direct consequence of Graf's addition formula (1.2). 

Acknowledgement. We thank the referee, A. P. Magnus and W. Van Assche for con- 
structive remarks concerning §4. 

2. Proof of the product formula 

In this section we present an analytic proof of the product formula (1.5). The proof 
uses two known formulas for the Hahn-Exton (/-Bessel function previously obtained by 
Koornwinder and Swarttouw [13] and Swarttouw [14]. 

The proof starts with the following formula, valid for \sxy\ < 1 and m G Z; 

m (s^xy' 1 ; g)oo(g m+1 ; g)oo ( q m s~ x xy- x , s~ x yx~ x _ 
(q m s- 1 xy- 1 ;q) 00 (q;q) 00 2 X \ q m+1 

— V" n z,,m+z (v 2 i 9)oo f . m +z+i\ z (x ; q)oo (®. nn z-\ 

^fr^ (9; 9)00 vr / (9; 9)00 

cf. [13, (4.5), (2.3)]. 

In (2.1) we take x = q2^ x ~ u + l \ y = q2i x+1 ) and s = q2 u+n for n G Z_|_ to get 

n w (r n ;i)x(f +1 ;t)x M {q m — n ,q-\ n nX+n+1 

q (q m — n ;qU(q;q)oo ^ { q m+1 ^ 

(2.2) 

= q z{n+1+l2V) Jx(q h(m+z) -,q)Jx-4<i l2Z -,q) 

z= — oo 

valid for m G Z, z/ G C, 9?(x) > —1 and all n G Z + . In (2.2) we use the notation (1.3). 
Multiply both sides of (2.2) by 

(q;q)oo (q u+1 ;q)n(q;q) n 

and sum from n = to 00. After interchanging the summations over z G Z and n G Z + , 
which is justified for 9fJ(x) > —1 and |_R| 2 < i ; c f. proposition A.l, we find 

q z Uq Um+z) ;q)J*-Aq^ z ; q)J v (Rq^ y+u+z) ,q) 

z= — oo 

^2 3) _ g^a: {q m+1 'i q)oo{q u+1 , q)oo i u ( u + y ) RU 

(9; 9)00(3; 3)00 



X 



\q ,_qhc { jj q 2 a n ^+^R 2n , - t q > q ■ a a x+n+1 ) 

l n m-v-n. n \ („v+l.„\ ( n -n\ * ^2^11 +1 ,</,</ I 
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We use for m G Z, n G Z + , 

,„_ m+ i.^ 



(g^" 1 " 1 ; 9)oo(9 ^ n ;9)oo _ / Am(m-l)-m(;y+7T,) (^ 



to see that the right hand side of (2.3) can be rewritten as 



( l) m g~^ m ^ '11)00(1 j 9)oo g\(v-m){v-rn+y)+\m{x+y) j^y 

(I] l)oo(i; l)oo 
( 2 - 4 ) 00 f_i)n |n(n+l) 



X 

n=0 



-.m — v—n „—n 



Now we apply the product formula for the Hahn-Exton 9-Bessel function, cf. [14, (3.1)], 
which is valid for a, 6, x, fx, v G C provided afoc 7^ 0. Explicitly, 



T / \ 7 / t \ (? ]l)oo{l^ ]l)oo v-in v + u, 

J v [ax;q)J tlL (bx\q) = r— r a \Tx +At 

(?; 9)00(9; 9)00 

(2 ' 5) ^ (-l)"(6x) 2n 9^^ +1 ) (q- n ,q- n ~» 

From (2.5) we see that (2.4) equals 

{-l)™ q -l™J m (Rq^+y)-q)J v _ m (RqUy+«-™ ). q ^ 

which proves the product formula (1.5). 

3. Proof of the addition formula 
The proof of the addition formula (1.4) uses the orthogonality relations 

00 

(3.1) i m+z Ui Uz+m) ;i)Ui ¥l+m) ;i) = s z ,i, 

m= — 00 

for zjEZ, > -1, cf. [13, (2.11)]. In [13, §3] it is shown that (3.1) can be viewed as 
a 9-analogue of the Hankel transform (or the Fourier-Bessel integral). 

Multiply both sides of the product formula (1.5) by q m J x (q^^ +rn ^] 9) for / G Z and sum 
over m from —00 to 00. We interchange the summations over m G Z and zGZ, which is 
allowed for U(x) > —1, |_R| 2 < 1; c f. proposition A. 2. An application of the 

orthogonality relations (3.1) and replacing m by — m yields 

j v (Ri^ v+1+v) ;i)J x -Ai ¥ ;i) 

00 

= (-i) m 9- im ^(^9^ (x+y) ;9)^ +m ( J R9^ (y+m+iy) ;9)^(9^°- m) ;9). 

m= — 00 

Since J- n (z;q) = (-l) n q? n J n (zq? n ;q), n G Z, cf. [13, (2.6)], the addition formula (1.4) 



:„ ,.„j 
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4. The limit case q j 1 of the addition formula 

In this section we present a (formal) limit transition of the (/-analogue of the addition 
formula (1.4) for the Hahn-Exton (/-Bessel function to Graf's addition formula (1.2) for 
the Bessel function. 

First we recall the (/-gamma function, cf. [3, §1.10], 

(4.1) r , W = J^(i- 9 )-, Bmr, W = r W . 

The g-gamma function can be used to see that (formally) 

(4.2) ]miJ v (z(l-q)/2;q) = J v (z). 

See [13, App. A] for a rigorous limit result of this type. 

In order to obtain Graf's addition formula (1.2) from (1.4) as q j 1 we have to consider 
the quotient of two Hahn-Exton (/-Bessel functions, which we rewrite as an infinite sum of 
quotients of two Hahn-Exton (/-Bessel functions of equal order. Explicitly, 

, 4 3 x Jx(q^ z ~ k) ;q) = q iu( z -k) (gj g)m m( i+i {x - v) ) J x -Aq^ z ~ k+m) ;q) 

J x _ v {q^ z ;q) ^ (q;q) m J x - U (q^;q) 

Note that the Hahn-Exton (/-Bessel functions on the right hand side of (4.3) only differ 
in the argument by a (half-) integer power of q. Equation (4.3) holds for J x _ v (q2 z ; q) ^ 
and 3?(x — v) > — 1. It can be proved by substituting the series representation (1.3) for 
the Hahn-Exton (/-Bessel function J x -v in the nominator on the right hand side of (4.3), 
interchanging summations and using the (/-binomial formula. 

The limit q j 1 in (4.3) can be taken with the help of the following proposition. The 
proof of the proposition is an easy adaptation of the proof of [17, thm. 1], which we will 
not give. 

Proposition 4.1. Suppose {p k (x;n) \ k G Z + ,n G N} is a series of functions satisfying 
a recurrence relation 

(4.4) x 2 p k (x;n) = a k+1 , n p k+1 (x; n) + b kj7l p k (x] n) + a fc)n p fc _i(x; n) 

with a k n > 0, b kjn G R. and initial conditions po(x;n) — f(x;n), pi(x;n) — g(x;n). 
Assume that the zeros of p k (x; n) are real for all k G Z + and all n G N and that 



(4.5) 



p fc _i(x;n) 



Pk(x;n) 



C 

< — , V/c,nGN, 
o 



uniformly for x G K for any compact subset K of C\R with d{K, R) > 5. Moreover, 
assume that 

lim a n n = A > 0, lim b n n = B G R, 

n^oo n^oo 
lim a k,n ~ a k-l,n = lim b k,n ~ &fc-l,n = 0, 
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Then we have 

(4.6) lim Pn + i^n) =p fx 2 ~B 



n^oo p n (x;n) \ 2 A 

uniformly on compact subsets of C\R with p(x) = x + \J x 2 — 1 (and the square root is 
defined by \p(x)\ > 1 for x [— 1, 1]). 

Next we apply proposition 4.1 to 
(4.7) Pk(x;n) = {-l) k q-^ k J 2na+ p{xq~^ k ;q) 

with q replaced by c x l n for a fixed c G (0, 1). We assume a > and (3 > —1, so that all 
zeros of pk{x;n) are real by [10, cor. 3.2]. It follows from the second order g-difference 
equation for the Hahn-Exton g-Bessel function, cf. [15, (17), (18)], that (4.4) is satisfied 
with 

a k , n = hkn = q\l + q^+?) 

and q = c 1 / 71 '. Then the conditions on the coefficients ak, n and bk, n of proposition 4.1 are 
easily verified. Moreover, A = c 1+a > and B = c(l + c 2a ) G R. It remains to prove that 
(4.5) holds. 

Lemma 4.2. With pk(x;n) defined by (4.7) for a > and f3 > —1, the estimate (4.5) 
holds. 

Proof. Consider the Wall polynomials for < b < 1, 



^ ;6;9) = ( -W^; 2V?1 (" b° ;q > x ) 



From [17, (2.5), cor. 1] we get for p G N 



w p - 



-i(x-Mq) < gVgj ~ 9 P )(1 ~ fegP" 1 ) < 9 p v^ 



w p (x;6;q) 

for all x G if', K' compact in C\[0, 1], 5' = d(K', [0, 1]). Replace b, x by q» +1 , y 2 q p - m and 
take the limit p — > oo. The Wall polynomials tend to the Hahn-Exton g-Bessel function 
uniformly, cf. [13, prop. A.l], and we get 



(4.8) 



Jjyq 2 m ;g) 



for all y E K, K compact set of C\R with 8 = d(K, R). The required estimate (4.5) for 
Pk{x;n) defined by (4.7) follows from (4.8). □ 

Consequently, by proposition 4.1, 

/„ ON ,• J2n a +f3(xC-^C-^,C^) 1 C(l + C 2a )-X 2 

(4.9) hm — ^ — — =p = p(a;;o!,c), /£?+- = r— • 
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In (4.3) we replace x, z, q by 2na, 2n7, c», and we take n — > oo. Iterating (4.9) and 
interchanging limit and summation shows that (4.3) tends to 



(4.10) c vl K -^c ma p k - m = c"V(l " ^p" 1 )- 



with p+l/p = (l + c 2a — c 27 )/c a , v < 1 and c 7 G C\R. Interchanging limit and summation 
can be justified using the estimate (4.8), [11, lemma A.l] and dominated convergence. 

Finally, divide both sides of the addition formula (1.4) by J x - V {q^ z ; q). Replace x, z, q 
by 2na, 2n7, c~ as before and y, R by 2ra?, R(l — q)/2. Take n — > oo and use (4.2), (4.10) 
to see that (1.4) formally tends to 



oo 

5] p k J k (Rc a+ ^J l/+k (Rc^) = c"^(l - p^cT^c^ 7 ), 

fc= — oo 



which is easily rewritten as Graf's addition formula (1.2) for the Bessel function. 



5. The limit case q j 1 of the product formula 

A (formal) limit transition of the product formula (1.5) for the Hahn-Exton (/-Bessel 
function to product formula (1.6) for the Bessel function is presented in this section. 

We first rewrite the product formula (1.5). Use J m (z;q) = (—l) Tn q~2 rn J- rn (zq~2 rn ;q), 
replace m by — m, and use the formula 

J x _„(z; q) = z~» f; { -J^^ HlHx) Q) 



k=0 



for > —1, and the notation J °° f{z)dm q {z) = J2T=-oo Q z f(<l^ z ) to write (1.5) as 
(5 1) 00 a\ r°° 

J2 W >Vk q Hl+$x) / z - VJ ^ zRq ^y+u). q) j x{zq -lrn. q )J x ( zq hk. q)dmq{z) . 

k=o yQ'iQh 



In order to calculate the limit of the integral on the right hand side of (5.1) we assume 
x = 2na : q = c~, for a > 0, c G (0, 1), so that we can use (4.7) (with (3 = 0) to find for 
r G Z + 

/•OO 

/ z 2r J x {zq-^ m - 1 q)J x {zq^ k -q)dm q {z) = 
(5-2) J ° 

(_i r +fc c £(m-fc) c -r c -£r / Z ^p n+m+k „)p„( Z ; n)rfm i (z), 

Jo c ™ 
by shifting the summation parameter. The orthogonality relations (3.1) imply 

/>oo 

/ p k (z: n)pi(z; n)dm i (z) = <5 fc /. 
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Repeating the first part of the proof of theorem 2 of Van Assche and Koornwinder [17] 
results in 



(5.3) z p w (z; n) Pfc+n (,; n)dm^ (*) = - y^ ^ _ = _ y dz 

for r,fc£ Z_|_, A = c 1+a , S = c(l + c 2q ) (as in §4) and T fc (cos#) = cos k9 is the Chebyshev 
polynomial of the first kind of degree k. 

In (5.1) we replace R, x, y and q by ^R(l — q), 2na, 2nr\, c«, and we use (4.2), (5.3), 
(5.2) to see that (5.1) (formally) tends to 



J u+m {Rc^)J m {Rc a+ v) = 

(5.4) 

k=0 



V (- u )k r ak ( ixfe+m [ B+2A (Z)- 1 ^ T (Rrn , T k+m {{z - B) /2A) 
h *U { } L-2A [ c } y~c } V 4A* -(z- BY 



as n — > oo. In the integral on the right hand side of (5.4) we replace (2 — B)/2A by cosV>, 
so that the right hand side of (5.4) equals 



(5.5) 27r 



- y { —^c ak (-ir +k r (2c a cos^ + 1 + c 2a )-^ 

fc=o fc! 

x J v (Rc 71 ^/2c a cos xj} + 1 + c 2a ) e "^ +m ^ d^- 
Replace -0 by -0 + 7r in (5.5) and use the binomial theorem to see that (5.5) equals 



— / J v (Rc v Jl + c 2a - 2c a cos ^ — '- —e-™* 

2tt J v v y (l + c 2 «-2c«cos^)^ 



(5.6) — / J v (Re Vl + c 2a -2c« cos V) — — e" <m ^ #. 



Equating the left hand side of (5.4) and (5.6) yields a formula equivalent to the product 
formula (1.6). 

Appendix. Justifications 

In the appendix we investigate the absolute convergence of two double sums occuring in 
sections 2 and 3. The estimates used are based on the estimates given in [13]. In particular, 
cf. [13, (2.4)], 

(A.l) \Mz;q)\<\z»\ [ -« ' qlZ Uq) °° 



(9; q)°o 

For integer order the Hahn-Exton g-Bessel function satisfies the estimate 
(A.2) \J n (z;q)\<\z\ ^- q >; q ^ q) ~ \\ ( „ 

by [13, (2.4), (2.6)]. A similar estimate also holds for the Hahn-Exton g-Bessel function 
with integer or half-integer power of q as argument because of the symmetry, cf. [13, (2.3)], 

(A.3) Ja(q^;q) = Mq^;q). 
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Proposition A.l. The double sum 

^ (?- +1 ;9)„(9;?),. ^ 5 1(4 ,5) ,? 



n=0 



is absolutely convergent for m G Z, > — 1, |_R| 2 q> 1 + 5f ( x )+ K (y) < 1. 
Proof. Assume m > and use (A. 2) and (A. 3) to estimate 

oo oo 

2 = 2 = 

for > —1, Vn G Z_|_ and C is a constant independent of n of which the values may 
differ with each occurence. Consequently, the part ^^L ^2T=o * s absolutely convergent 
for M(.r) > -1. 

Similarly we estimate 

— m oo 

(A.4) ^ \(f^ 1+ * v) J x (q* {m+z) \q)Jx- v (q* z \q)\ < C ^ ^^"-^M)^^-!). 

2= — oo 2=m 

Shift the summation parameter to run from to oo, estimate q z ( z ~^ by q^ z( - z ~^ /(q; q) z 
and use the o</?o summation formula to see that (A.4) can be estimated by 

Qq-nm ^_q$l(x-is)-n+l+m. ^ q q —^n{n—l) q nU{x—u) 

So the part Y^=q Y^=-oo * s absolutely convergent for \R\ 2 q 1+ ^( x )+^(y) < 1. 

The remaining sum over z is finite and causes no problems. The case m < is essentially 
the same. □ 

Proposition A. 2. The double sum 

oo oo 

q m Uq w+m) ;q) E q z Uq h(m+z) -,q)Jx-Aq l2Z ;q)MRq^ y+u+z) ;q) 



is a 



bsolutely convergent for dt(x) > —1, \R\ 2 q 1 +^( x ')+^(y') < 1. 



Proof. First we estimate the inner sum over z as a function of m. For the Hahn-Exton 
g-Bessel function of order x and x — v we use the same estimates as before, and from (A.l) 
we get 



\J v (Rq* (v+u+z) ;q)\ < Cq^ 



1, z > 0; 

\ R \-2z q -zmy+v) + l) q -±z(z-l) ^ z < 



where C is independent of m. Denote by S m (l) the inner sum over z from max(0, — m) to 
oo, then we find, using these estimates, for 9ft(x) > —1, 



\S m (l)\ < C\ _ , 



g imK(x) 5 m > 0; 
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For S m (2), the inner sum over z from — oo up to min(0, — m), we obtain 

{^m(iSR(a;)+SR:(y))|^|2m^ m > 0; 
? -i m »(x)+|m(m-l) ) m < 0. 

The remaining finite sum over z from min(0, — m) + l to max(0, — m) — 1, denoted by S m (3), 
can be estimated by 



\S m (3)\<C< 



^-imK^)^ m<0. 



To estimate the double sum of the proposition we combine (A. 2) and (A. 3) to obtain 
an estimate on J x {q^^ l+m ^]q). The sum m = — oo to min(0, — /) and the finite sum from 
m = min(0, — /) + 1 to max(0, — /) — 1 are absolutely convergent. The sum m = max(0, — /) 
to oo is absolutely convergent for $l(x) > — 1 and q 1 +^( x )+^(v) |_R| 2 < i. □ 
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